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Abstract 

The inclusive gluon jet function is evaluated at two-loop accuracy. This function is 
relevant for resummations of large perturbative logarithms in collider processes involving 
low-mass gluon jets. The jet function corresponds to the imaginary part of the gluon 
propagator in light-cone gauge, which is adopted for the calculation. In addition to the 
leading jet function, the power-suppressed two-gluon jet functions are given and their 
renormalization is discussed. 



1 The gluon jet function 



Cross sections for the production of low-mass jets at large momentum transfer often factorize 
into hard, jet and soft functions [HI2]- While the hard functions encode the physics associated 
with the large momentum transfer Q, the jet functions describe the propagation of collinear 
partons inside the jets. The soft functions mediate low energy interactions between jets and 
have a low characteristic scale M 2 /Q, where M is the typical invariant mass of the jets. 

Soft-Collinear Effective Theory (SCET) implements this factorization on the operator level 
[3J IH E] . In this framework the hard functions appear as Wilson coefficients of operators built 
out of soft and collinear fields. At leading power, the operators do not contain soft fields and 
the soft interactions can be decoupled from the effective Lagrangian. After the decoupling, 
the soft functions are given by matrix elements of light-like Wilson lines along the directions 
of large energy flow. The matrix elements of the collinear fields associated with the final state 
define the jet functions. The collinear fields associated with the initial state at hadron colliders 
define (generalized) parton distribution functions, also called beam functions [5J. By solving 
the renormalization group (RG) equations of the hard, jet and soft functions, one can resum 
contributions to jet cross sections which are logarithmically enhanced by ratios of the different 
scales. 

Most of the early applications of SCET were concerned with processes mediated by electro- 
weak currents with only two directions of large energy flow. Such processes involve only quark 
jets. The corresponding jet function was computed to two-loop order already some time ago 
[7]. The gluon jet function first appeared in computations of radiative T decays [HI [9], [10] , but 
more generally, it appears in any process involving gluon jets. This includes any jet process at 
hadron colliders and the production of three or more jets at e + e~ machines. The jet function 
is given by the vacuum matrix element of two gluon fields, 

d 4 xe^ <0| Al{x)A b M |0> = Y,W^\P ~ Px) <0| ^(0) \X)(X\ A b M\0) . (1) 

x 

Loosely speaking, it measures the probability that a gluon field produces a jet of particles 
with momentum p from the vacuum. The field A^x) is the gluon field in light-cone gauge. It 
is related to the field A^x) in a general gauge via 

A"{x) = A a ^{x)t a = W\x) [iD»W(x)} . (2) 

The light-like [n 2 = 0) Wilson lines 



W(x) = P exp ^ig s J ds n ■ A(x + sn)\ (3) 

ensure that the jet function is invariant under gauge transformations which vanish for x — > oo. 
The symbol P indicates path ordering, and the conjugate Wilson line W* is defined with 
the opposite ordering prescription. In SCET the jet function is defined as a matrix element 
of collinear fields which describe the propagation of energetic particles with small invariant 
masses p 2 <C (n ■ p) 2 . However, after decoupling the soft fields, the collinear Lagrangian is 
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equivalent to the QCD Lagrangian. For the computation of the matrix element (|T|) we may 
therefore use QCD Feynman rules. 

The one- loop result for the gluon jet function was given in [11], where it was used to perform 
soft-gluon resummation for direct photon production at next-to-next-to-leading logarithmic 
(NNLL) accuracy Without too much effort it should be feasible to push the logarithmic 
accuracy for vector boson production processes even one order higher. Based on the results 
of [121 H3] (see also [1^1 [15]). the necessary three- loop anomalous dimensions were derived in 
The two- loop hard functions are also known and can be extracted from the results of 
[T6| ITT] . The only missing ingredients to perform a N 3 LL resummation for W, Z and photon 
production at large transverse momentum are the two-loop results for the soft function and 
the gluon jet function. In the present paper, we compute the gluon jet function at two-loop 
order. 

In addition to the inclusive jet function considered in this work, it is also interesting to 
study jet functions that incorporate kinematical restrictions on the final state X which appears 
in ([T]). Such jet functions become relevant when a jet-algorithm is applied to the final state. 
The one-loop quark jet function for Sterman- Weinberg jets has recently been computed in 
|18j . and both quark and gluon jet functions for general recombination and cone algorithms 
have been worked out in [19] (for related work cf. [201 EI])- 

The most general parameterization of the gluon jet function consistent with the identity 
n ■ A{x) = and invariance under a rescaling of the light-cone vector reads 

J d 4 xe^(0|^(xK(0)|0) 

= 5 ab g 2 s 9(p°) 

where the strong coupling constant g s on the right-hand side is the bare coupling. In the 
construction of SCET one introduces a conjugate light-cone vector n M with n ■ n = 2 and 
decomposes the gluon field as 

A»(x) = At(x) + — n ■ A(x) , (5) 

with n ■ A±(x) = n ■ A±(x) = 0. The jet function J g (p 2 ) then corresponds to the matrix 
element of two transverse fields Aj_(x). Since the component n ■ A(x) is suppressed in the 
limit p 2 <C (n-p) 2 , only the function J g {p 2 ) is relevant in leading-power factorization theorems. 
However, independence of the reference vector ensures that the field components always 
appear in the combination ([5]) and for completeness we will compute the full two-point function, 
including the power-suppressed piece K g {p 2 ). 

In the following section, we will discuss the technical aspects of the two-loop computation 
and present the bare, unrenormalized, two-loop results. As a check of our calculation, we have 
performed the computation both in Feynman and light-cone gauge. The renormalization of 
the jet functions is discussed in Section [3j As an interesting aside, we discuss the choice of 
the operator basis necessary to renormalize the power-suppressed jet function K g {p 2 ). 



-9^u + 



n ■ p 



Up 2 ) 



(n ■ pY 



K 9 {P 2 



, (4) 
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2 Two-loop calculation 



In order to calculate the jet function, it is convenient to rewrite (jl]) as the imaginary part of 
the time-ordered product 



d*xe^{0\T{A;(x)A b u (0)} \0) = S* £ 



-9»u + — — Jg{P ) + ,_ _\v lCg(P 



n ■ p 



(n ■ p) 



with 



Up 2 ) = -im [ij g (p 2 )] 

7T 



K g (p 2 ) = -Im[i)C a (p 2 )] . 

7T 



(6) 
(7) 



This form shows that the jet function is given by the imaginary part of the gluon propagator 



in light-cone gauge n ■ A{x) = 0, where the Wilson lines become trivial and A h 



x 



g s A^(x). 



For dimensional reasons, the function J g (p 2 ) ~ (p 2 )~\ while K, g {p 2 ) ~ (p 2 ) . Thus only J g (p 2 ) 
receives contributions from single-particle intermediate states, while JC g (p 2 ) starts at one-loop 
and has a cut for p 2 > but no pole. 

Since the Wilson lines are absent in light-cone gauge, it is convenient to use this gauge in 
the computation of the jet function. The free gluon propagator in light-cone gauge is given by 



p 2 + iO 



n ■ p 



(8) 



Note that we are not adopting the Mandelstam-Leibbrandt (ML) prescription to regulate the 
n ■ p — » singularity. The ML prescription 
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^ (9) 

n ■ p n ■ pn ■ p + zO 

cures the collinear singularity in the propagator, but in SCET this singularity has a phys- 
ical meaning. The Wilson line W(x) and the associated light-cone propagators arise from 
expanding QCD diagrams around the large-energy limit and the choice of the zO-prescription 
is dictated by the QCD diagrams. The loop integrals contributing to the jet function are un- 
ambiguously defined in dimensional regularization. They depend on a single four momentum 
p. If they involve k light-cone propagators they scale as (n ■ p)~ k . Since the dependence on 
n -p is analytic, the final result of the integral is independent of the sign of the zO-prescription 
adopted for the light-cone propagator. 

As a check of our result we performed the calculation both in light-cone and Feynman 
gauge. Both gauges offer some advantages, but overall the light-cone gauge computation is 
more efficient. The number of diagrams in this gauge is significantly smaller, because no 
Wilson line diagrams and no ghost-loop contributions need to be included. At one-loop the 
light-cone gauge calculation reduces for instance to the first two diagrams in Figure [TJ At 
two-loop order, the light-cone gauge calculation involves the diagrams shown in Figure El 
while the number of diagrams is again more than doubled in Feynman gauge. Furthermore, 
in light-cone gauge, standard tools for QCD computations can be used, e.g. we generate the 
diagrams with FeynArts [24J. The only disadvantage of light-cone gauge is the proliferation of 



3 



Figure 1: Diagrams contributing to the gluon jet function at one- loop. In light-cone gauge the 
ghost contribution, represented by the third diagram, as well as the diagrams in the second 
row, which involve emissions from the Wilson lines (denoted by the crosses), are absent. 



light-cone denominators. In the diagrams in Figure [2J each gluon propagator has an associated 
light-cone denominator, while in Feynman gauge such denominators only arise from Wilson 
lines. 

The most general two-loop integral that appears in the light-cone gauge calculation has 
the form 



Jdk Jdl ^ ^ _ ^ + ^ ^ + ^ ^ + l + p ^ b3 

1 



(n ■ k) Al (n ■ l) A ' 2 (n ■ k — n ■ l) Ai {n ■ k + n ■ p) Bl {n ■ I + n ■ p) B ' 2 (n ■ k + n ■ I + n ■ p) Bi 

(10) 

This class of integrals contains more master integrals than those which appeared in the two- 
loop computation of the quark jet function [7]. Nevertheless, after using partial- fractioning 
identities and symmetry relations, no additional master integrals remain after solving the 
integration-by-parts relations [25] for the integrals which actually appear in the diagrams. To 
solve the integration-by-parts relations we use Laporta's [26] algorithm as implemented in the 
public codes AIR [27] and FIRE [2B] as well as an independent code written by one of us. 

After combining the irreducible two-loop diagrams with one-particle reducible ones, we 
find that both calculations yield the same bare result. We obtain 



^ 2 > = Z^{l + f^(z^ l Me) 

rrl 1 / 1 \ 2e 



< (47r)2 ( -pf-io ) K J ^( £ ) + C A n f T F J Af (e) + C F n f T F J Ff + n 2 f T F J ff (e)] j , (11) 

where a s refers to the MS coupling constant, which is related to the bare coupling constant 
a° s via Z a a s /j, 2e = e- e7B (4vr) £ a° with Z a = 1 - /3 a s /(4vre) and /3 = 11/3 C A - 4/3 Tpn/. The 
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Figure 2: Two- loop diagrams that arise in the light-cone gauge calculation. 



one-loop coefficient is [TT] 

Me) ... e 4 
and the two-loop coefficients are found to be 



n f T F 



^r(2-6) 2 r(6) 

T(4 - 2e) 



8 55 1 / 7T 2 152 \ 1 / 184C 3 llvr 2 3638 



3e 3 



23tt 4 1496C 3 161tt 2 57415 



6 



27 



180 



9 



20 188 1 /2tt 



3e 3 9e 2 



+ 



27 162 ' 

536 \ 400C 3 , 74tt 2 12880 



9 



9 



27 



2 _ 55 

- + I6C3 - T 
e 3 



J //( e ) = 77^ 



16 160 8tt 2 
9^ + 27^-^7- + 16 



The result for the bare function JC g takes the form 



Ki{e) + 



2c 



(12) 



(13) 



(An) \-p 2 - lO J w (4tt) 2 \-p 2 - tO 

Only a single color-structure appears at one-loop 

4T( 2 - e)r(-e)r(e) 



[Cl^(e) + C A n/T F i<: A/ (e)] 



-C 4 e e7£ - 



T(2 - 2e) 



(14) 
(15) 



and only two structures at two-loop level 

K AA (e 



19 148 1 / 151C 3 7tt 2 1082 



e 4 ' e 3 



3e 2 



6 



9 



17vr 4 428Cs 41vr 2 7672 
~72 3 9" + ~W 



K Af (e) 



4 38 298 44C 3 13tt 2 219 



3e 2 9e 



9 



27 



(16) 



The bare jet functions Jg e (p ) and K^ are (p 2 ) then follow by taking the imaginary part using 



-im u-p 2 - *ori = -0(p 2 )^M(p 2 )° 

7T 7T 



(17) 



The function Jg are (p 2 ) is a distribution in p 2 whose explicit form is obtained after expanding 



1 f p 



p 2 



1 00 

V)+E 



n=0 



rv. 



p2 



18) 



The star-distributions are generalizations of plus-distributions to dimensionful variables, their 
definition can be found in [29J. 



3 Renormalization 

3.1 The leading jet function J g (p 2 ) 

In momentum space the leading jet function renormalizes via the convolution 

. 2 

J g (p 2 , fi) = f dp' 2 Z Jg (p 2 - p' 2 , p.) J^(p' 2 ) . 
Jo 



(19) 



Both the jet function and the Z-factor are distribution valued, see (Tl8|) . For this reason it is 
more convenient to perform the renormalization in Laplace space, where it is multiplicative. 
The Laplace transformed function j g is defined as 



■ f 1 Q2 



dp 2 e- up2 J g (p 2 ,{i)i with 



Q 2 e^ 



(20) 



The function j g is also what is needed to perform soft-gluon resummation in the momentum 
space formalism of [30]. It fulfils the RG equation 



di^( ln f'^ 



Q 2 
/i 2 



-2T A In — - 2 7 J « 

cusp 2 ' 



jg (In^-,// 



(21) 



where T^ nsp denotes the cusp anomalous dimension, which is known at the three-loop level [31] . 
The jet anomalous dimension r y Jg was inferred to three loops in [TT] using RG invariance of 
the direct photon-production cross section and the three-loop results for the hard anomalous 
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dimensions 



the quark jet anomalous dimension [52], and the Casimir scaling property 



of the soft function. Expanding the anomalous dimensions as r^ sp = $^L o r„ V4T 



Y^=Qln (f^) n+1 , the two-loop solution to the RG equation f[2"TT) takes the form 



a. 



■ «'£)■ 



«) 2 y + (-ft + 3o^)r4 



+ (if + (7o') 2 " /W* + cfT?) y + (T! J ' + lt'c{- - l3 c{')L + el 



(22) 



The coefficients entering at this order are 

'67 7T 2 



1 



7o' 



4C A , Tf = 4Cj 



9 



20 



9 



T F n f 



(23) 



-A), Ti =C. 



, , 1096 llvr 2 
1 -— + — + 160 



27 9 

Together with the one-loop coefficient [TTj 

'67 2tt 2 



27 



9 



C 



.4 



9 



20. 
~9 



TpTlf 



(24) 



the solution (I2"2"j) completely determines the two-loop jet function up to the constant c 2 9 , which 
we compute in this work. 

In Laplace space the jet function renormalizes multiplicatively, j g = Zj g jg are , and Zj g fulfils 
the same RG equation (121]) as the renormalized jet function. Solving this equation [T2l [13], 
one derives the following expression for the logarithm of the Z-factor: 



09 4^ 



F A 1 



-^-7(^1"^- ■ 



4e 3 27 2 



Q 2 



7^ ( ^In^ + 7o J9 ) - % + i ( Tfln^ + 1? 



Q 2 



(25) 



Since all the necessary anomalous dimensions are known, the Z-factor is completely deter- 
mined at the two-loop level. The cancellation of all divergences l/e n , for n = 1...4, in 
the renormalized result provides a strong check of our calculation. We finally obtain for the 
non-logarithmic two-loop coefficient 



C 2 —C\ 



( 20215 362vr 2 88C 3 17tt 4 



V 162 



27 



36 



1520 1347T 2 I6C3 
+ C A n f T F I — — + 



27 



27 



(26) 



55 



+ C F n f T F [ ~ + I6C3 ) + n)Tl 



400 8tt 2 
~8T ~ ~27 
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Numerically, for nj = 5 flavors, this yields 

J g (L, n) = 1 + (|^) (-2.961 - 7.667 L + QL 2 ) 

+ (-pj (—58.58 + 44.39 L + 82.46 L 2 — 61.33 L 3 + 18 L 4 ) , (27) 
which may be compared to the quark case 

] q (L, //) = 1 + (^) (0.560 -4L + 2.667L 2 ) 

+ {^j 2 (—36.34 + 32.14 L + 43.25 L 2 

Not surprisingly, the corrections are larger in the gluon case. On the other hand, for a s = 0.1 
they are numerically small. For \L\ < 1, the two-loop corrections amount to ±0.4%. To gauge 
their phenomenological relevance, we have included the corrections into the resummed result 
for the direct photon cross section [11]. At the Tevatron, and for = 50 GeV, the change in 
the cross section due to the inclusion of the two-loop gluon jet function is around +1%, while 
the two-loop corrections from the hard and soft functions (which are, however, only partially 
known) are both around +5%. 

3.2 The subleading jet function K g (p 2 ) 

While there exists only one physical operator at leading power, we need to include additional 
power-suppressed operators to perform the renormalization of the subleading jet function. 
This is obvious, since K g {p 2 ) vanishes at tree level, but is divergent at one-loop order. To 
discuss renormalization, it is convenient to work with the scalar operators 

Oj{x) = ^T{Al{x)A h M) , 

(29) 

O k {x) = T{8- A a {x) d ■ A b (0) } + DOj(x), 

whose vacuum matrix elements give rise to the jet functions J g (p 2 ) and K g (p 2 ), respectively. 
In d — 4 dimensions, Oj(x) ~ l/x 2 and 0/c(x) ~ l/^ 4 , so the second operator is power- 
suppressed with respect to the first one. 

For our purposes it is sufficient to concentrate on power-suppressed two-gluon operators, 
which mix into operators with more than two fields, but not vice versa. Starting from the lead- 
ing operator Oj, infinitely many power-suppressed two-gluon operators can be constructed. 
Examples include 

1 (x) = -DO J (x), 2 (x) = 2 \ - Oj(x) , 0,{x)=A^Oj{x). (30) 

— X z + l\J ox 



— 17.48 L 3 + 3.56 L 4 ) . (28) 
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In Laplace space one finds that the vacuum matrix elements of these operators are 



Ji(Q 2 ) = e^ Q* — 3g (Q*) = 0(a s ), 

UQ 2 ) = e» (Q 2 Y ! Q &Q l2 {Q' 2 Yl 9 {Q' 2 ) 
Jo 

UQ 2 ) = # B Q 2 ] 9 (Q 2 ) = e>* Q 2 + 0(a s ) , 



e ,E Q 2 
1 - e 



+ 0(a s ) , 



(31) 



where j g is the leading (bare) jet function. Notice that, at tree-level, j\ vanishes and the 
difference j'2 — J3 also vanishes in d — 4 dimensions. Since the operators Oi, O2 and O3 
have the same field content, including spin structure, their physical matrix elements cannot 
be distinguished and it is therefore sufficient to include an arbitrary linear combination with 
non-vanishing tree-level matrix element in the operator basis. If several of such operators are 
included, the renormalization scheme can always be chosen such that only one combination 
has a no n- vanishing renormalized matrix element. In this sense the additional operators 
are similar to evanescent operators which appear in matching computations onto four-quark 
operators [331 . However, in contrast to this case, we are not forced to include the additional 
operators to obtain a basis which closes under renormalization. 

From fl3Tl) it is obvious, that the simplest choice corresponds to only including the operator 
O3. In momentum space, the corresponding jet function is just the integral over the leading 
jet function, 



For completeness, let us also give the momentum space representation of 2 , which reads 



Note that the integral in the second term is not ultraviolet convergent in d = 4 dimensions 
and can thus not be performed on the level of the renormalized jet function J g (p 2 ,fi). A 
similar behavior has been observed in convolution integrals relevant for power corrections to 
the process B — > X s 7 [34J. Subleading quark jet functions have recently been studied in 
|35j . where it was noticed that an operator basis which includes only the equivalent of 0\ is 
insufficient to perform renormalization. It was further speculated whether the operator basis 
closes under renormalization if O3 is added. Our discussion makes it clear that this is the 
case, as far as additional operators related to the leading jet function are concerned. 

Let us now discuss the renormalization of the two subleading jet functions k g and j 3 . The 
renormalized jet functions are obtained as 




(32) 




(33) 




(34) 
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The first row follows from the fact that j'3 = e lE Q 2 j g and that the leading operator j g cannot 
mix into the subleading operator. Of the two other coefficients only Z%j enters at one-loop 
order, and we find 

Z k j = , (35) 

71 e 

and 

UQ\ = V* (l - In OP) . (36) 

To obtain the renormalized result at two loops, we would need to compute the renormalization 
factor Zkk at one-loop order. To separate out corrections associated with Zkk from the ones 
proportional to Zjy, one would need to compute a matrix element which is sensitive to the 
spin of the gluon fields in the operators. The simplest possibility is to compute a one-gluon 
matrix element of the operator 0% at one-loop. 

4 Conclusions 

We computed the gluon jet function at two-loop order. Performing the calculation in both 
Feynman and light-cone gauge, we argued that the latter leads to considerable simplifications. 
The current calculation yields one of the two missing ingredients to perform the N 3 LL soft- 
gluon resummation for W, Z and photon production at large transverse momentum. In 
addition, we computed the subleading two-gluon jet functions and discussed how to construct 
a minimal basis of power-suppressed operators which closes under renormalization. 
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